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Abstract : We consider an electron coupled to a random distribution of 
point vortices in the plane (magnetic impurities). We analyze the effect of 
the magnetic impurities on the density of states of the test particle, when 
the magnetic impurities have a spatial probability distribution governed by 
Bose or Fermi statistic at a given temperature. Comparison is made with the 
Poisson distribution, showing that the zero temperature Fermi distribution 
corresponds to less disorder. A phase diagram describing isolated impurities 
versus Landau level oscillations is proposed. 
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1. Introduction 

The problem of a 2 dimensional electron gas coupled to a static random magnetic field 
has been a subject of interest in the past few years [1]. Particular attention has been 
paid to localisation properties of such systems. In the case of Gaussian disorder with zero 
mean, all states seem to be localised [2]. On the contrary, they are delocalised in the case 
of an uniform magnetic field. Therefore, the question arises about the role played by a 
mean-field description of a random magnetic field [3,4]. 

The model we are interested in consists of a planar electron of electric charge e coupled 
to a random magnetic field 9 . What we mean by random magnetic field is a distribution 
of infinitely thin vortices carrying a flux (ft, modelizing some sort of magnetic impurities, 
characterized by the dimensionless Aharonov-Bohm coupling a = e(ft/2n (i.e. (ft in unit of 
the quantum of flux). This system is periodic in a with period 1 and since there is no 
privileged orientation of the plane, it is invariant by changing a into —a, implying that a 
can be restricted to the interval [0, 1/2]. 

In order to study the effect of statistics on the disordered magnetic impurity systems, 
one may evaluate perturbatively in a the average one electron partition function, i.e the 
Laplace transform of the average density of states. 

In previous works [4], we focused on magnetic impurities obeying a Poisson distri- 
bution (which actually corresponds to the Bose case at zero temperature, or at infinite 
temperature). We observed in particular a transition for a c ~ 0.35 between an almost free 
density of states for an isolated impurities system (a > 0.35) and an oscillating Landau 
like density of states (a < 0.35). 

In this letter, we will consider the random magnetic impurities system as a gas of 
particles of a given density p, with a distribution obeying Fermi or Bose statistics at a 
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temperature T v . We will first properly define the perturbative expansion of the average 
partition function. Then, we will explicitly compute at order a 2 contributions to the av- 
erage partition function. We will argue that in the Fermi case, at zero temperature, the 
average density of states always displays Landau like oscillations, implying that there is 
no transition between an isolated impurity disordered phase and a mean magnetic field 
phase. This result is not unexpected, since Fermi statistics clearly tends to homogeneize 
the impurity configurations, leading to a less disordered situation. 

2. The Model 

2.1. General Formalism 

Let us consider an electron coupled to a random magnetic field given by a distribution 
p(r) of magnetic impurities. This means that p(r)dr is the number of impurities at position 
r in the infinitesimal volume d 2 r. The Hamiltonian is given by 

where we have explicitely taken into account the coupling of the magnetic field to the spin- 
up (+) or down (— ) degree of freedom of the electron endowed with a magnetic moment 
f 1 = ~ 2m (thus an electron with a gyromagnetic factor g=2). 

In the case of a discrete distribution p(r) = J2i^( r ~ r i)> where the index i indices 
the impurities, the spin-term is a sum of contact terms. It corresponds to a choice of a 
peculiar self adjoint extension [5]: in the (+) case, the wave functions vanish at the location 
of the impurities (hard-core boundary condition), whereas in the (— ) case singular wave 
functions are considered at the location of the impurities (attractive-core conditions). In 
order to extract the short distance behaviour of the wave functions, a non-unitary wave 
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function redefinition has been used [6]: 

N 

^(r) = nir-r i | ± >±(r) (2) 

i=i 

The generalisation of this transformation to a continuous distribution p(r) is 

^± (r) = e ±^Vp(Oin|r-r^± (r) (3) 

The Hamiltonian acting on ^(r) rewrites 

H+ = --d z d- z - — / dz'dz'^-^-d z (4) 
m m J z — z 

H~ = --d z d- z + — [ dz'dz'^-^-d, (5) 
m m J z — z' 

where the complex coordinates in the plane have been used z = x + iy, d z = \{d x — id y ) 

and dzdz = d 2 r. H or H can be used indifferently to compute the partition function, 

since it is by definition the trace of a function of H. In the sequel we will concentrate on 

the spin up coupling, keeping in mind that the spin down analysis could be easily done 

following the same lines. 

Up to now p(r) has not been specified. If a Poissonian distribution is choosen, p(r) is 

defined by its cumulents 

p(ri)...p(r fe ) = pS(ri - r 2 )5(r 2 - r 3 )...(5(r fc _i - r fc ) (6) 

Here however, we deal with quantum statistics for the impurities themselves, so p(r) has 
to be defined as 

p(r)=V + (r)V(r) (7) 

with 

V'(r) = ±f d 2 ka(k)e-^ r (8) 
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y,+ (r) = — / d 2 ka+(k)e ikr (9) 

a + (k) and a(k) are the creation and annihilation Fock space particle operators , with the 
commutation rules 

[a(k),a+(k')] = <5(k-k') (10) 

for bosons, and 

{a(k),a+(k')} = <5(k-k') (11) 

for fermions. 

Thus, the impurities, considered as a quantum gas, have a temperature T v and a chemical 
potential p, which determines their mean density p. The average over disorder of an 
operator Q consists in 

<Q>= 3 Ji l (12) 

where the impurity second quantized Hamiltonian 

H v = J d 2 k{^ ~ M)a+(k)a(k) (13) 

describes the equilibrium of the impurity gas in the grand-canonical ensemble. Note that 
we consider here quenched impurities, which are in thermodynamical equilibrium. Note 
also that the Poissonian distribution dP(ri) = dri/V can be seen as the particular case of 
Bose distribution at T v = 0. The impurities indeed condensate in the zero energy iV-body 
wave function V>( r i> r 2j • • • ,*n) = (-j=) N , leading to the iV-impurity Poissonian distribu- 
tion dP = V>*(ri,r 2 , . . • .rviriri.ro. . . . , r N )dr 1 dr 2 . ..dr N . 

We wish to evaluate perturbatively the average one electron partition function at 
inverse temperature (5 

,2a s „ [P .„ fPp- 1 

p=0 



< r|e _/3if |r 



^ la fl 3 rPp-i 
>= E(— y / d &- / df3 p G^ l3l {v,v l )i, + {v' 1 )^{v\ 
~n fn Jo Jo 
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r) B 



. _ r -G /3l _ /32 (r 1 ,r 2 )..^ + (r' p )^(r' p ) F ^-G^(r p ,r) (14) 

where integrations over the position variables r, and r'j are implicit. G/3(r, r') is the free 
electron propagator 

c " lri ' r! ' = 5 e " 5 ""'" 1 ' < i5) 

Averaging over disorder yields expressions like 
which can be evaluated using the contractions 

d 2 k 



±gf(r,r')=<4;(r)i> + (r')>= J ^(1 ± (16) 
g±(r,r') =< ^(r)^(r') >= J ^n^~^ (17) 



rik stands for the Bose-Einstein (upper sign) or Fermi-Dirac (lower sign) distributions. 

nk = e /3„(fcV2m-M) T 1 

One has the relation 

gf(ry)=g±(r,r')±6(r-r>) (19) 

To summarize, the perturbative expansion of the one electron average partition func- 
tion can be represented in terms of Feynmann diagrams given by rules quite analogous to 
those of finite temperature second quantised formalism: 
electron line: Gp{ri,Tj) 
forward impurity line: 9f( r i> r 'a(i)) 
backward impurity line: 9^{ Y 'ii v 'a{i)) 

impurity loop: gf (r-, r •) = p 

electron-impurity vertex: ^ x , d Zi 



For a given diagram of order p, the electron propagates from its initial to its final position r 
via r p , r p _i,...ri, the location of electron-impurity interaction, with temperatures 0, f3 p ,..., 
(3\, [3. The p vortices located at position r^... r' p , undergo a permutation a. At each a(i), 
it corresponds a vortex line: 

backward line if a(i) > i 

forward line if a(i) < i 

and a loop if a(i) = i. 

In the Fermi case, each diagram is affected by the signature of the permutation a. 
The dimensionless parameters at work are the rescaled average density X 2 p in unit of the 
electron thermal wavelength A 2 = 2iT(5/m, the rescaled average density X 2 p in unit of the 
impurity thermal wavelength A 2 = 2wP v /m, and the Aharonov-Bohm coupling constant 
a. 

Clearly, one expects that in the limit (3 V — > 0, i.e. Boltzmann statistics, with uncor- 
rected randomly dropped impurities, one recovers Poisson distribution. Also, as already 
advocated above, one expects that in the limit (3 V — > oo, in the Bose case, one has again the 
Poisson distribution, whereas the Fermi distribution leads to a less disordered situation. 
In the sequel, one will concentrate on the relative interplay between the dimensionless 
parameters X 2 p, A 2 /), a to study the phase diagram of the magnetic impurity system. 

2.2. Mean-field expansion 

Consider first diagrams that are entirely built by impurity loops (Fig. la). These di- 
agrams do not involve the many-body statistical correlations of the impurity distribution, 
and are thus independant of the statistic. Therefore, they yield the same contribution as 
in the Poisson case [4] . 
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The order a p (p > 0) term writes 

Summation over p yields as it should the partition function per unit volume of the mean 
magnetic field e < B >= 2npa (i.e. in the mean-field limit a — > 0,p —* oo, pa finite) 

Z<b> = l—^-pl^i) (20) 

2tt 2sinh/3^^ yy h 2m ' y ' 

The global positive shift in the Landau spectrum is a direct manifestation of the hard-core 

boundary conditions on the wavefunctions 

2.3. Second order expansion 

Non trivial diagrams (Fig. lb), i.e. not mean-field diagrams, appear at second order in 
a. Let us denote the fugacity by z = =pe^" M . In the interval z g] — 1, 1] one has the 
expansion 

OO 2 

n(k) = ± {-z) n e~ n ^h (21) 
n=l 

As a result 

oo 

g f (Tj) = S(t - r') + £(-*)»G ft/ j B (r,r') (22) 

?1=1 

OO 

5fe (r,r')=±E(-^) nG, ^( r ' r/ ) ( 23 ) 

n=l 

The density p is related to z by 

p = ^ln(l + z) (24) 
Using (||,|D, the dia gram of Fig. lb has the contribution 



_ _ ™? g / ,> is "\ (25 ) 

2 4/3, ^ n + p I 7o x„p + x(l-xW 



At high temperature (T v > T), in the Boltzmann limit, i.e. z —* both in the Bose and 
Fermi cases, one gets D(z) — > pa 2 /2, which is precisely the Poisson distribution result. 
The correction reads 



D{z) = ^ r (l- , -\lp) 
At low temperature (T v < T), D(z) rewrites 



D(z) = EEL 



■^%-^ ( .2 xg (O^i^^m, h q - m (z)h m+ i(z) 



(26) 



(27) 



where the functions h q (z)'s are defined in the interval ] — 1, 1] by 



h q (z) = J2 



i-zr 



n=l 



(28) 



The low temperature limit corresponds to z — > +oo in the Fermi case and z — > — 1 in the 
Bose case. It happens that (p8|) can be analytically continuated in ] — 1, +oo[ by noticing 
that 

ho(z) = --4- (29) 



and by using the recursive relation 



l + Z 



h q (x) 
X 



dx 



(30) 



As a result, D(z) in (27) is valid in the interval z E] — 1, +oo[. It follows that in the Fermi 
case and in the low temperature limit T v — > 0, one can use the expansion 



(In*)* 9 E ^ 2) , rn (M^ 



n=l 



-l) n 1 



(31) 



which is valid for z > 1. 
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Which contributions does D{z) yield in the low T v limit? In the Bose case, some care 
is required, since in ( |27|) the limit z — > cannot be interchanged with summations. We 
checked numerically that ( f27[) indeed yields pa 2 /2, i.e. the Poissonian result as expected. 
One obtains for the average partition function Z at order a 2 

Z = 1(1 - \x 2 pa + ^(A 2 p)V(^ + ^) + • . .) (32) 

In the Fermi case, using the identity YXX C m=o (^r^yy (^TT)T = (q+i)f( g !)^ ■ (H© 
yields 

lim D(z) = ^![l--L(l + A 2 p-e- A2 ' , -2yAV/' / ^dye-^) (33) 

At this point one can consider either a low impurity density limit X 2 p « 1, or a high 
impurity density limit X 2 p » 1. At low density X 2 p << 1, one finds that the average 
partition function per unit volume rewrites as 

Z = A^ (1 " \ X2pa + \ {X2p)2a2{ W P + h X2p) + ' • (34) 

At order a 2 , the p 2 a 2 term is missing, a situation quite different from the Poissonian case 
(j3^), where this mean-field term precisely dominates in the mean-field limit. Therefore, in 
the Fermi case at low temperature, the low density expansion is not adapted to describe 
the mean-field limit. On the other hand, at high density X 2 p 3> 1, (1331) leads to 



z = > - + k*>wk + aWS - w P + + • " <35) 

where the leading mean magnetic field term is indeed present. 

3. Fermion case at zero Temperature: the Landau regime 

3.1. The ordered phase 
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Before the system reaches the mean-field limit (pO[), we expect H an intermediate 
regime caracterised by smooth Landau oscillations in the spectrum. This intermediate 
regime is identified as the ordered phase by opposition to the phase with no oscillation 
(disordered phase). The order a 2 is quite instructive to give information on the way the 
system reaches the mean-field limit. Consider indeed the case of Poissonian impurities (|3^). 
We observe 1/p corrections to the mean-field at any order of the perturbative expansion in 
(ap) n . On the other hand, ( |35| ) shows corrections to the mean- field of order 1 / (p^/p). This 
implies that the system approaches more rapidly its mean-field limit when the impurities 
are fermions at zero temperature, rather than Poissonian. In other words, the system 
is less disordered, since a Fermi distribution of impurities is more homogeneous than a 
Poissonian one. 

Let us generalise these considerations at any order a n of perturbative theory. One has 
to evaluate < p(ri)p(r2)...p{r n ) >, which can be rewritten as 

n -. „ n 

< p(n)..p(r n ) >= E E ~< / dr' 1 ..dr' p p(r' 1 ,..,r' p )l[5(r q -r' m ) (36) 

is the set of all possible surjections from (l,..,n) to (l,..,p) and p(ri,..,r p ) is the 
p-body correlation function 

p(r 1 ,..,r p )= E e(^)5ft(ri - r a{1) )...g b (r p - r a(p) ) (37) 
In the case of Fermions at zero temperature, one has the correlator 



1 



g b (v) = -J^J l (^F P r) (38) 
r V 7T 



For example, using 



< p(ri)p(r 2 ) >= p 2 - ^ qjxtv^lri - r 2 |)] 2 + P 8{v x - r 2 ) (39) 

n - r 2 r 7r 
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altogether with (14), yields the contibution (33), in addition to the mean- field term. 

At high density (X 2 p 3> 1), because of Fermi exclusion, the p-body correlation function 
becomes 

p(n, .., r„) = f?[l - - £ 5(r t - rj ) + 0( J-)] (40) 
P i<j p yP 



When (40) is used for evaluating < p(r\) p(r2) ...p{r n ) >, one indeed finds that corrections 



to the mean-field term (pa) n are of order l/(p^/ r p). 

3.2. Absence of a pure disordered phase 

We have just seen that corrections to the average magnetic field limit are less important 
in the Fermi case. Could it be that the statistics of the impurities alter the occurence of 
the transition itself [Q]? First let us remark that the average partition function can be 
expressed as 

Z = ± I F(\ 2 p,\ 2 v p,a) (41) 

which means that the average density of states is a function of E/p, X 2 p and a. This is 
due to the fact that <#,(r) is p times a function of ^fpr and \\p. In ( |14| ) together with 



(pq) , rescaling /3j into /3j//3, rj into r,/A and into ^/pr^, immediately leads to (41). In 
particular at T v = 0, X 2 Z is a function of X 2 p and a. Considering the expansion 

X 2 Z = 1 + ^a(a - l)X 2 p + c 2 (a)(X 2 p) 2 + c 3 (a)(X 2 p) 3 + ... (42) 

we can show that C2 = 0. Use simply 

and conclude in general that p{r±, ..,r p ) is at least of order pP. In particular 

/ 9(r 1 ,r 2 )= /3 2 -[ 5;) (r 1 -r 2 )] 2 (44) 
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starts at order /? 3 , whereas p(r) = p. Therefore, no contribution of order p 2 is to be 
found in < p(ri)p(r2)...p(r n ) >, which means that the average partition function does not 
contain terms of order p 2 at any order a n . 
The specific heat H 

C = k P 2 -$2 l * Z ( 45 ) 

gives a lower bound for the critical value a c at which the system does not exhibit any 
Landau oscillations. More precisely, if the correction to C — Cq (Co = k) is negative when 
(3 — ► 0, then oscillations are already present since, at small /3, 

C - C = 2Af r r dEdE> i< ^ V>d< *fJ/ V> (E - E-f + . . . (46) 
jo dE dE' 

Here, the small (5 expansion ( ^4|) yields 

C - Co = -Ifc(^) 2 a 2 (l - a) 2 + . . . (47) 
4 m 

which is always negative. Therefore the average density of states always displays Landau 
oscillations. 

4. Conclusion 

To complete this analysis, let us emphasize again that for intermediate magnetic im- 
purity temperature, the average partition function has been shown in ( j4l| ) to scale as 1/A 2 
times a function of X 2 p, X 2 p and a. This implies that the average density of states is a 
function < p(E/p, A 2 /3, a) >. Since Poisson distribution is recovered in the Boltzman limit 
T v — > oo (since then p(r±, ..,r p ) — * pP), one interpolates between the Poisson and the zero 
temperature Fermi cases simply by varying X 2 p from to oo. When A 2 /? = 0, there is a 
transition at a c ~ 0.35, whereas when A 2 p = oo, we have just shown that no transition 
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occurs at all. Therefore, we expect that, for A^p sufficiently small, a transition still occurs 
at a critical value a' c > 0.35, and for X^p sufficiently big, no transition occurs anymore, 
meaning that the system is always Landau like, in the whole interval a € [0, 1/2] (Fig 2). 
It is not clear if the transition observed in the Poissonian case can be interpreted as a 
phase transition, possibly related to a localisation-delocalisation transition. But magnetic 
impurity distributions do influence this transition by actually reordering the system when 
the density correlations are increased. 
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